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Introduction 



The dynamically cold component of the Kuiper Belt is host to a population of very widely separated, 
near-equal mass binary systems. Such binaries, representing the tail of the separation distribution of the 
more common, more tightly-bound systems, are known to have on-sky separations up to ~ 4". Their 
wide separations make them highly valuable due to their delicacy and sensitivity to perturbation, and 
also makes them relatively easy targets to characterize from the ground. Parker et al. (2011) present a 
ground-based characterization of seven such systems with separations at discovery ranging from 0."5 — 
4", and we will adopt these systems as the prototypes for the ultra-wide binaries of the Kuiper Belt. 
Here we present the prospects for using future large-scale ground-based optical surveys (with LSST 
as our baseline survey) to measure the orbital properties of a large sample of these widely separated 
Trans-Neptunian Binaries (TNBs). 

1. Seeing distribution for LSST observations 

The seeing at the LSST site is well measured, and the optical characteristics of the LSST imager 
have been extensively modeled. A summary of the estimated delivered seeing is presented in the LSST 
Science Book (LSST Science Collaboration 2009), and Table 1 contains the FWHM of the delivered 
PSF in the r and i bands. We adopt these percentiles for determining the repeated resolvability of TNB 
systems. 
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a continuous function PDF for the see- 
ing distribution in order to generate synthetic observational circumstances. We adopt a PDF for the 
measured seeing distribution at the LSST site which closely reproduces the peak and 25th, 50th, and 
75th percentiles of the r and i band (see Table 1). 

2. Projected binary separation 



The fraction of observations which are resolved by LSST will drop for binaries with smaller primary 
radii, given the same Hill sphere separation distribution as a function of radius. This is because the Hill 
spheres of TNOs scale linearly with their radii, so a TNB with a smaller primary radius will have a 
smaller on-sky separation for a physical separation of a given fraction of its Hill sphere. The observed 
fraction of binaries as a function of primary radius must be corrected for the fraction of a Hill sphere 
which is resolvable as a function of primary radius in its discovery survey. Such characterization is 
critical, as the trend of binary fraction with primary radius is a strong indicator of the collisional history 
of the Kuiper Belt (Petit & Mousis 2004, Nesvorny et al. 2011, Parker & Kavelaars 2012), providing an 
orthogonal constraint from size distribution measurements. 

Assuming that the binary system is distant and on a circular heliocentric orbit, so r ~ r®, the 
projected Hill radius can be estimated for the following two cases: given the primary and secondary radii 



Rp and R s , when R P > R s , the projected Hill radius R H ~ 18".3(p /l g cm 3 ) 5 100km), while in the 
case that R P ~ /? s , /?# ~ 23".l(p/l g cm 3 ) 5 100km). 

The second case, where the secondary is of comparable mass to the primary, will be used throughout 
the rest of this paper. 

Grundy et al. (2008) found that the most probable projected separation of TNB systems at random 
epochs (for random orientations and a variety of eccentricity distributions) was within one percent of the 
mutual semi-major axis (a m ), with 84% of the observed separations being within a factor of two of a m . 
Therefore for our analytical estimate we adopt a m as the characteristic separation from which to estimate 
the on-sky separation. 



Table 2 demonstrates the implica- 



Table 2: Radii of resolvable equal-mass TNBs 



0.05 56 km 67 km 77 km 

0.10 28 km 33 km 36 km 

0.20 14 km 16 km 19 km 



tions of the variation of Hill radius size a ^ R " Rp for which system is resolved in. 

changing with primary radius, illustrat- _ 25th percentile 50th percentile 75th percentile 
ing the primary radius at which a binary 
with a given a m jRu must be in order for 
that separation to translate into an on- 
sky separation sufficient to be resolved under 25th, 50th, and 75th percentile conditions at LSST For 
smaller primary radii, only objects occupying larger fractions of their Hill spheres will be resolved by 
the survey. 

At present six known TNBs fall in the 0.05 < a m /Rn < 0. 1 range (Parker et al. 201 1 , Shepherd et al. 
2011), while just three binaries fall in the 0.1 < a m /Rji < 0.2 range, and only one TNB is known with 
a-m/RH > 0.2 (the exceptional system 2001 QW322). 



2.1 Separation of known ultra-wide TNBs 

In addition to these analytical estimates, we 
simulate the observations of known wide TNBs in 
order to determine the fraction of LSST observa- 
tions which will resolve systems with compara- 
ble properties. The binaries studied in this work 
are drawn from Parker et al. (2011), which rep- 
resent the widest binaries with characterized or- 
bits. These seven systems were found to occupy 
between 7% and 22% of their mutual Hill spheres, 
with physical separations of several tens to one 
hundred thousand kilometers. Mutual eccentric- 
ities varied from ~ 0.2 to 0.9. 

Given the baseline of LSST observations of 
ten years, we sample each system's mutual orbit 
over ten year segments of that system's heliocen- 
tric orbit centered at random epochs. Over these 
ten year segments, we make synthetic "observa- 
tions" of each system by determining the observed 
separation of the system at 100 random epochs, 
given the Earth-binary separation and the angu- 
lar orientation of the system. This process is re- 
peated for hundreds of random ten- year segments 
for each binary system. 

At each synthetic observation, the LSST seeing distribution was used to generate an estimate of the 
seeing at the time. If the binary separation at the epoch exceeded the FWHM of the sampled seeing, the 
binary was said to be resolved during this epoch. Figure 1 illustrates the results of this exercise. We find 
that all seven of the mutual orbits in Parker et al. (201 1) would be resolved in some fraction of observa- 
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Figure 1: Fraction of LSST observations for which 

binary systems with mutual orbits drawn from 
Parker et al. (2011) are resolved — eg., only ~ 5% 
of all binaries with properties comparable to those 
in Parker et al. (201 1) are resolved in fewer than 
10% of visited epochs by LSST — the remaining 
95% of wide TNBs will be sampled in more than 
10% of the total number of observations. 



tions from the LSST site, only 5% of all possible orientations of these seven systems will be resolved in 
fewer than 10% of all images taken from the LSST site — with the other 95% of possible orientations 
resolvable in more than 10% of the same images. Given the baseline of 230 r-band observations per field 
over the 10 year LSST survey, this would net > 23 resolved epochs for 95% of these ultra-wide TNBs. 

3. Estimates of sample size 

The estimated H-magnitude completeness for LSST TNO observations is greater than 90% com- 
pleteness (differential) for H < 7.5 (LSST Science Collaborations 2009). The CFEPS survey (Petit et al. 
2011) estimate that in the cold classical Kuiper Belt there are approximately 4,000 objects brighter than 
H ~ 8, and that the luminosity function of this population is steep with a ~ — 1.2. This would indicate 
roughly 1,000 objects brighter than H ~ 7.5 which would be recovered in the highly complete LSST 
sample. Extending to the 50% (differential) completeness magnitude for LSST, H ~ 8.6, we would have 
to extrapolate the results of the CFEPS survey over more than half a magnitude. Adopting the same 
luminosity function slope over this range, we estimate approximately 20,000 TNOs brighter than this 
magnitude; because of the steepness of the luminosity function, most detections will be in the last few 
luminosity bins and we can conservatively approximate the differential completeness of 50% as the cu- 
mulative completeness and adopt an estimate of 10,000 TNO detections from the cold classical Kuiper 
Belt brighter than H ~ 8.6. 

The brighter, more complete sample also conveniently corresponds to a radius of ~ 77 km (when 
adopting a geometric albedo of 10%). This is the radius for which even binaries occupying just 5% 
of their mutual Hill sphere will be resolvable in the 75th percentile conditions of the LSST site (Table 
2). Given the literature estimate of a wide binary fraction of ~ 1.5% (Kern et al. 2006, Lin et al. 
2010), we would expect roughly 15 of the 1,000 TNOs detected by LSST in this bright sub-sample to 
be resolved binaries. Such binaries, being bright and easily separated, will have high signal-to-noise 
measurements of their mutual astrometry in many visits, and consequently their mutual orbits will be 
very well characterized by the LSST observations alone. 

The fainter, less complete sample will return far more resolved binaries. The limit of H ~ 8.6 corre- 
sponds to a radius of approximately 46 km. This radius is slightly smaller than the limit resolvable in the 
25th percentile seeing at LSST given an a m /Ru < 0.05. However, this is larger than the limit resolvable 
in even the 75th percentile seeing for a m jRu ~ 0.1. As a crude estimate, four of the 10 known wide 
binaries have separations in excess of 0.1/?//, so we can estimate the fraction of binaries in this separa- 
tion class to be 0.4 x 0.015 ~ 0.006. We then expect 60 resolved binaries in this class found by LSST, 
given our estimate of 10,000 cold Classical TNOs detected in total. Such binaries will have relatively 
low signal to noise observations as they are near the photometric limit of LSST, but they will be resolved 
under most conditions and orbital characterization will be possible. 

If we consider the limiting primary radius which can be resolved in 25th percentile conditions by 
LSST given a separation in excess of 0.05/?//, we find Rp ~ 56km which translates to H ~ 8.2 given an 
albedo of 10%. This magnitude limit is comparable to the magnitude at which the CFEPS estimates are 
defined, with approximately 4,000 brighter objects in the cold classical Kuiper Belt in total. Given a wide 
binary fraction of 1.5%, we expect 60 resolved binaries in this sample. While these objects will have 
moderate signal to noise observations, they will be resolved under a smaller number of circumstances 
(~ 25%) than the previous two cases. 

The preceding sub-samples have some overlap, and when accounting for the variation of observable 
fraction of Hill sphere with radius along with the increasing number of TNOs in the sample with de- 
creasing radius, the sum total of TNBs expected to be resolved in at least the 25th percentile conditions 
(and therefore with > 60 resolved observations in r band alone) is approximately 100. This binary sam- 
ple represents the systems that will be orbitally characterized by LSST itself; given the rapid increase 
of binary fraction with decreasing mutual separation (eg., Kern & Elliot 2006) and that for eccentric 
systems the observed separation can significantly exceed the projected mutual semi-major axis in some 
fraction of observations (eg., Grundy et al. 2008), a significantly larger number of systems will likely be 
discovered by LSST but not directly orbitally characterized. 



4. Value of Large Mutual Orbit Sample 



• Exceptionally detailed dynamical decomposition: Do prograde TNBs exhibit a marked difference 
in mutual orbit properties compared to retrograde? Such a difference might indicate separate 
formation pathways dominating for prograde and retrograde orientations (Schlichting & Sari 2008, 
Nesvorny et al. 2010, Parker et al. 2011). 

• Detailed host population statistics: Does the Kernel component of the cold Classicals (Petit et al. 
2011) or the collisional family of Haumea (Brown et al. 2007) host wide binaries? Are there 
any such wide binaries amongst the hot Classicals? Is the one resonant wide equal-mass binary 
(Sheppard et al. 201 1) an anomaly, or does it represent a population of resonant wide TNBs? 

• Is the fraction of binaries correlated with primary radius, indicating collisional modification of the 
Kuiper Belt (Nesvorny et al. 2011, Parker & Kavelaars 2012)? 

The final point is particularly important, as a conelation of binary fraction with radius would be 
particularly interesting from a dynamical standpoint; however, it would also significantly undercut the 
probability of detecting a large sample of wide binaries with LSST Sheppard et al. (2011) present a 
survey which discovered 1,000 faint TNOs and only find one system (the brightest in their sample) to 
be a binary. If such a low wide-binary fraction is characteristic of the small objects in the cold Classical 
population, then the likely sample size of wide binaries discovered by LSST will be drastically reduced; 
the estimates of the number of binaries in the bright subset discussed previously would be the least 
affected. However, the binary 2000 CF105 has a primary radius of ~ 32 km, making it extremely sensitive 
to collisional disruption. Its continued existence argues against a strongly collisionally evolved Kuiper 
Belt (Parker & Kavelaars 2012) and given that its albedo is high (~ 0.2 — 0.3) relative to the mean of the 
cold Classicals it may represent a bright member of a large population of small Kuiper Belt binaries. 

Conclusion 

The population of ultra-wide TNBs are valuable tracers of the dynamical history of the outer So- 
lar System, but at present the conclusions that can be drawn from them are limited by small-number 
statistics. Barring any strong correlation of binary fraction with primary radius, LSST will characterize 
roughly 100 ultra-wide TNBs and discover many more. Such a sample will allow substantially stronger 
conclusions to be made regarding the origins of these systems and their implications for the history of the 
outer Solar System. Careful debiasing will allow disentangling of albedo distribution and binary fraction 
trends, further constraining the collisional grinding of the Kuiper Belt. 
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